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Abstract 

Learning- and memory-related processes are thought to result from dynamic interactions in large-scale brain networks that 
include lateral and mesial structures of the temporal lobes. We investigate the impact of incidental and intentional learning 
of verbal episodic material on functional brain networks that we derive from scalp-EEG recorded continuously from 33 
subjects during a neuropsychological test schedule. Analyzing the networks' global statistical properties we observe that 
intentional but not incidental learning leads to a significantly increased clustering coefficient, and the average shortest path 
length remains unaffected. Moreover, network modifications correlate with subsequent recall performance: the more 
pronounced the modifications of the clustering coefficient, the higher the recall performance. Our findings provide novel 
insights into the relationship between topological aspects of functional brain networks and higher cognitive functions. 
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Introduction 

During the last years evidence has accumulated suggesting that 
an improved characterization of time-variant interactions between 
different regions within the complex network brain can be 
achieved with graph-theoretical approaches (see [1—6] for an 
overview). Within this framework a network (or graph) is 
considered as a set of nodes (or vertices) and a set of links (or 
edges) connecting the nodes. Functional brain networks can be 
derived from direct or indirect measurements of neural activity 
(e.g., electroencephalogram (EEG), magnetoencephalogram 
(MEG), or functional magnetic resonance imaging (fMRI) data). 
The sampled brain regions are usually considered as network 
nodes, and network links represent interactions between pairs of 
brain regions that can be assessed by evaluating interdependencies 
(see [7-1 1] for an overview) between their neural activities. The 
resulting connection schemes can then be characterized by 
network metrics [12] such as the average shortest path length or 
the clustering coefficient. The average shortest path length L is the 
mean number of steps along the shortest paths for all possible pairs 
of nodes. The clustering coefficient C is the mean of the local 
clustering coefficients of all nodes and quantifies the tendency of 
nodes to form local clusters. The local clustering coefficient of a 
node is the fraction of triangles among all connected triples with 
the node as their center. Large values of both L and C are 
characteristic for an ordered, lattice-like structure; low values of L 



and C are typically observed for random networks. When 
analyzing functional brain networks, the clustering coefficient is 
often interpreted as a measure of the local efficiency of information 
transfer, and/or of the robustness of the network to deletion of 
individual nodes. Similarly, the average shortest path length is 
often interpreted as a measure of the global efficiency of a network 
to transfer information between nodes [2,5]. 

Higher cognitive functions are thought to result from dynamic 
interactions of distributed brain areas operating in large-scale 
networks [13-18]. Using graph-theoretical approaches, previous 
studies provided evidence that cognitive tasks-such as motor 
learning or (n-back) working memory-specifically modify func- 
tional brain networks derived from recordings of neural activities 
(EEG, MEG, fMRI) acquired during the tasks [19-35] and that 
these modifications appear to be associated with task performance 
[19,24,29,35]. Other studies reported on associations between 
properties of functional brain networks recorded during a resting 
state condition and cognitive performance metrics assessed 
independendy during neuropsychological evaluation [36-39] . 

Most of the aforementioned studies investigated cognitive 
functions that mainly involve prefrontal areas [40] and primary 
association cortices but only rarely temporolateral neocortex as 
well as important structures of the mediotemporal lobe, namely 
hippocampus and surrounding perirhinal cortex and parahippo- 
campal gyri. The lateral and mesial temporal lobe structures are 
known to be involved in learning- and memory-related processes 
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[41-49], and their functionality has frequently been investigated 
with fMRI [50-53], intracranial electrophysiological techniques 
[54-59], or via non-invasively recorded event-related electric 
potentials or magnetic fields [60]. 

Here we investigate whether incidental and intentional learning 
of verbal episodic material differentially modifies functional brain 
networks and whether modifications are related to subsequent 
recall performance. To this end, we assess— in a time-resolved 
manner-global statistical characteristics (L and C] of functional 
brain networks that we derive from ongoing multichannel EEG 
data recorded non-invasively from healthy subjects and from 
epilepsy patients during a neuropsychological test schedule. 
Chronic epilepsy and its treatment is known to impair cognitive 
processes and to induce functional reorganization and behavioral 
compensation [61-63]. Investigating both healthy subjects and 
epilepsy patients thus provides a spectrum of recall performances 
that is required to study possible relationships between learning- 
induced network modifications and memory. 

Materials and Methods 

Ethics statement 

The study was approved by the local Ethical Committees, and 
all subjects gave written informed consent. 

Subjects 

Thirteen patients (age mean of 36 + 1 1 years; 8 females; mean 
duration of epilepsy 19 + 9 years) with pharmacoresistant epilep- 
sies of suspected temporal or extratemporal neocortical origin as 
well as twenty healthy controls (age mean of 32 + 8 years; 10 
females) were included in the study. All patients had been 
submitted for pre-surgical evaluation at the University of Bonn 
Epilepsy Program [64] . Seven patients were diagnosed as having a 
unilateral temporal lobe epilepsy (mesial seizure-onset zone in four 
patients and lateral seizure-onset zone in three patients). In five 
patients the seizure-onset zone was located in neocortical 
structures (central, parietal, parieto-central, parieto-occipital), 
and one patient was diagnosed as having multiple seizure-onset 
zones. Surgery was performed in eleven patients (selective 
amygdalohippocampectomy; extended lesionectomy) and led to 
complete seizure control in six patients. 

Neuropsychological test schedule and EEG recordings 

Before participating in the neuropsychological test schedule, 
subjects were informed about its general outline. They will be 
asked to learn a set of words (learning task) and later on, to 
remember the learned words (retrieval task). They will also be 
asked to perform an extra control task during which they just have 
to listen to words, however they were not informed about a 
subsequent retrieval task. The test schedule thus consisted of two 
blocks of tasks which involved incidental (block 1) or intentional 
(block 2) learning and retrieval of verbal material (cf. Figure 1): 

• Block 1: during a learning task (Tj) subjects were orally 
presented 15 highly frequent German nouns (cf. [65]) five 
times in randomized order. During the retrieval task (7T) 
subjects were asked to freely recall words they listened to 
during T[. 

• Block 2: during a learning task (T^) subjects were required to 
learn 15 highly frequent German nouns (different from those 
presented during 7"') that were orally presented five times in 
randomized order. During the retrieval task (7^) subjects were 
asked to freely recall words memorized during Tl. 



Learning and retrieval tasks lasted 3 minutes each and were 
separated by a figural fluency task (to minimize recency effects, to 
discourage verbal rehearsal, and to avoid interference of any 
verbal information with the learned verbal material) of 1 minute 
duration followed by a rest phase of 3 minutes duration. Blocks 
were separated by a rest phase (with eyes open) of 5 minutes 
duration, and a baseline recording (7),) of 5 minutes was 
performed with eyes open preceding the experiment. In order to 
minimize possible order effects, we balanced the sequence of the 
blocks over subjects. 

The number of correctly retrieved words (N\ and N2I excluding 
perseverations) during the retrieval tasks served as measure of 
verbal memory (recall performance). Control subjects recalled 
9.1 ±3.5 (range: 3-14) words (Aq) during T[, and 1 1.5±2.8 (range: 
6-15) words (Nj) during T[. Epilepsy patients recalled 5.6±3.6 
(range: 0-11) words (N\) during T[, and 7.7 ± 3.5 (range: 3-13) 
words (/V2) during T' 2 . The correlation between Aq and N2 
(Pearson's correlation coefficient; control subjects: 
p = 0.6,p<0.01; epilepsy patients: p = 0.63,/><0.05) as well as 
the difference between the distributions of performance measures 
(t-tests; both groups: /><0.05) proved significant. Control subjects 
performed significandy better than epilepsy patients (t-tests; both 
tasks: /><0.02), and both groups showed a significantly increased 
recall performance during T 2 (t-tests; both groups: /><0.05). 

During the whole examination procedure EEG data were 
acquired continuously at a sampling rate of 254.31 Hz (16 bit A/ 
D conversion) within a bandwidth of 0-50 Hz from 29 electrodes, 
and right mastoid served as physical reference. Locations and 
nomenclature of these electrodes are standardized by the 
American Electroencephalographic Society [66]. In addition, we 
recorded right horizontal and vertical electrooculogram. In order 
to minimize the influence of technical (e.g., amplifier resets) and 
physiologic artifacts (e.g., eye movements/blinks or head move- 
ments) we applied a wavelet-based correction scheme as described 
in detail elsewhere [67]. 

Patients were investigated, on average, seven days prior to pre- 
surgical evaluation with blood levels of the anticonvulsant 
medication within the therapeutic range. 

Construction of functional networks 

To construct functional networks from EEG data we associated 
network nodes with EEG electrodes (K = 29) and inferred network 
links by estimating interdependencies between EEG time series 
from pairs of brain regions, regardless of their anatomical 
connectivity. For this purpose, and motivated by an increasing 
evidence for (phase) synchronization between distant brain areas 
to play an important role in learning- and memory-related 
processes [49,68,69], we employed the concept of time-variant 
phase coherence [70] which is an established method for studying 
time-variant changes in phase synchronization. The mean phase 
coherence [71] 



N 



expi(® m (j)-<f>„(j)) 



(1) 



is the temporal average of the differences of the instantaneous 
phases O of time series from nodes m and n, and N denotes the 
number of data points. By definition, r mn is confined to the interval 
[0,1] where r mn = 1 indicates fully phase-synchronized systems. We 
used the analytic signal approach [72,73] with which instanta- 
neous phases are obtained from the Hilbert transform of a time 
series. An important property of the analytic signal approach is 
that the the instantaneous frequency-in case of two or more 
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Figure 1 . Neuropsychological test schedule. The schedule consists of two blocks of tasks which involve incidental (block 1 ) or intentional (block 
2) learning and retrieval of verbal material. Each block is preceded by a rest phase of 5 minutes duration. Learning tasks (T[ and T' 2 , 3 minutes 
duration) are followed by a figural fluency task (FF, 1 minutes duration) and a rest phase of 3 minutes duration. Retrieval tasks (T[ and T' v 3 minutes 
duration) consist of a free recall of words listened to/memorized during the respective learning task. 
doi:1 0.1 371 /journal.pone.0080273.g001 



superimposed oscillatory components— relates to the predominant 
frequency in the Fourier spectrum [74,75] which may be subject to 
fluctuations in the time series. In such a case the instantaneous 
frequency varies rhythmically around the predominant frequency 
resulting in spurious estimates of the instantaneous phase. 
However, such effects can be reduced by taking the temporal 
average (cf. Eq. 1) or by applying adaptive approximation methods 
proposed recently [76]. From an electrophysiological point of 
view, we consider it more reasonable to look adaptively (via the 
Hilbert transform) at synchronization between predominant 
rhythms in the EEG than to look at synchronization in some a 
priori fixed frequency bands (e.g., via wavelet) for which there is no 
power in the time series (cf. Refs. [75,77,78]). 

For further analyses and to allow an interval-based estimation of 
the mean phase coherence, we split the EEG time series into 
consecutive time intervals of 16.1 s duration each (corresponding 
to N = 4096 data points) and estimated, for each time interval, the 
elements r mn of the phase synchronization matrix R. From this 
matrix we constructed binary networks using a thresholding 
approach, and set the non-diagonal elements of the adjacency 
matrix A to a mn = 1 if the corresponding entry r mn of R exceeded a 
threshold 8, and to a mn = 0 otherwise (a mm = 0 Mm). The choice 
of a threshold is not trivial, and there is currently no commonly 
accepted method to infer links from estimators of signal 
interdependence. However, since the number of nodes and the 
link density can affect network characteristics [79-85] precautions 
must be taken to minimize such influences when comparing 
between different networks. In order to determine 8 we requested 
that the networks are sparse (k«K) and that they do not 
disintegrate into unconnected network components as this leads to 
infinite values for the average shortest path length. Connectedness 
is guaranteed for mean degrees k > In (K), at least for random 
networks [86] , and with the fixed number of nodes K = 29 we 
obtained a mean degree k = 4, such that In (K) <k«K holds. 

Computation of statistical network characteristics 

Here we considered the average shortest path length L and the 
clustering coefficient C as global statistical characteristics of a 
functional brain network. We used an algorithm proposed in Ref. 
[87] to determine the shortest paths between all pairs of nodes 
from the adjacency matrix A. The length of the shortest path L mn 
between two nodes m and n is the minimum number of links to 
traverse in order to get from m to n (or vice versa). With this 
definition we derived the average shortest path length L of the 
network by taking the average of all L mn (m # n) with L mn < go . 
The clustering coefficient C m for node m can be defined as (see, 
e.g., [1]): 



C m = ( k m {k m - 1) E./ a "" a >^» *™ > 1 (2) 
V 0, k m e{0,l} 

where k m denotes the degree, i.e., the number of in- or out-going 
connections of node m. We averaged C m over all nodes to obtain 
the clustering coefficient C of a network. 

In the following, we consider the baseline recording (T/,) and the 
learning tasks (T[ and T^j only. We omitted the retrieval tasks (T[ 
and T^) and the figural fluency task in order to avoid the influence 
of movement-related artifacts. With the interval-based estimation 
of interdependencies we obtained for each subject 1 1 intervals for 
T[ and for T l 2 and 18 intervals for Tj. Averaging over intervals 
resulted in mean network characteristics which we denote as L, 
and C t (the suffix t denotes task). 

Results 

Task-related modifications of functional brain networks 

In Figure 2 we show, for an exemplary epilepsy patient and an 
exemplary control subject, temporal evolutions of the network 
characteristics during different neuropsychological tasks. Indepen- 
dent of tasks, both characteristics exhibited considerable fluctua- 
tions over time, and the patient's fluctuations were more 
pronounced. In both subjects, the learning tasks {T[ and Tj) did 
not appear to modify the average shortest path length of the 
functional brain networks, as we could not observe clear-cut 
differences in L between any tasks. The same holds true for the 
clustering coefficient of the patient's functional brain network. For 
the control subject however, the intentional learning task (Tj) led 
to a significant increase of the clustering coefficient C (when 
compared to both Tt, and T\; t-tests; /><0.05). Both network 
characteristics attained slighdy higher values in the control subject 
during all tasks. When comparing the characteristics' distributions 
from the patient group and the control group however, we could 
not observe any significant differences between the groups despite 
their significantly different recall performances. For the binary 
networks considered here, this observation is in line with our 
previous results [67]. In the following, we present our findings 
obtained from merged groups. 

In order to investigate whether the aforementioned influences of 
neuropsychological tasks on network characteristics extend beyond 
exemplary data, we estimated-for all subjects and each task- 
statistical properties of the distributions of L and C and evaluated 
possible differences between tasks (ANOVA with Greenhouse- 
Geisser and Huynh-Feldt corrections for departure from spheric- 
ity; /><0.05). The distributions differed between tasks for the 
clustering coefficient but not for the average shortest path length. 
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Figure 2. Exemplary temporal evolutions of global statistical network characteristics during different neuropsychological tasks. 

Time courses of average shortest path length L (left) and clustering coefficient C (right) from an epilepsy patient (red lines) and a control subject 
(black lines). Baseline recording (Zj); incidental learning task (T[); intentional learning task (Tj). Lines are for eye-guidance only. 
doi:1 0.1 371 /journal.pone.0080273.g002 



From the box-plots shown in Figure 3 it can be deduced that both 
learning tasks did not significantly modify the average shortest 
path length L of the functional networks. The intentional learning 
task T 2 , however, led to a small but significant increase of the 
networks' clustering coefficient (C T i > Cp and C T i > Cr b \ t-tests; 
/><0.05). 

Relationship between network modifications and recall 
performance 

We investigated whether network modifications induced during 
the learning tasks (T[ and Tj) are related to the number of words 
(Ni and JV2) subjects recalled during the subsequent retrieval tasks 
(T[ and T£j. Given the high interindividual variability (cf. Figure 3) 
we focused on the relative deviation of the average shortest path 
length and of the clustering coefficient during the learning tasks 
from the respective values during the baseline recording Tb, i.e. 
L* = (L t — Lx b ) I Lx b and C* = (C, — Cj-J / Cj b , where t denotes 
7j and T 2 , respectively. 

For the whole group of subjects, both characteristics deviated up 
to 20% from the baseline (cf. Figure 4). On average, however, 
deviations of the average shortest path length were close to zero 
(0.3% for T\; 0.2% for T 2 ; no significant difference), while those of 



the clustering coefficient amounted to 2.9% for T 2 , being 
significantly higher than the deviation for T[ (0.2%; t-test; 
/><0.05). Interestingly, only for the intentional learning task-for 
which subjects were explicidy instructed to learn the presented 
material-could we observe a significant positive correlation 
(Pearson's correlation coefficient p = 0.40;/) = 0.03; cf. Figure 5) 
between the relative clustering coefficient C* , and the number N2 

of subsequendy retrieved words. Due to the aforementioned 
significant correlation between the recall performances, C* , was 

also positively correlated to the number of words N\ acquired 
during T[ (Pearson's correlation coefficient p = 0.37;/? = 0.05). Of 
note, there were no significant correlations between the recall 
performances and the characteristics L and C of the networks 
from the baseline recording T/,. 

Conclusions 

We investigated whether functional brain networks are differ- 
entially modified during incidental and intentional learning of 
verbal episodic material and whether modifications are related to 
subsequent recall performance. We assessed the two global 
statistical characteristics-clustering coefficient and average shortest 
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Figure 3. Comparison of the distributions of global statistical network characteristics from 33 subjects for different 
neuropsychological tasks. Box-plots of the average shortest path length L (left) and the clustering coefficient C (right) for the baseline 
recording (Tj) and for the incidental (T[) and the intentional learning task (Tj). Bottom and top of a box are the first and third quartiles, and the band 
inside a box is the median. The ends of the whiskers represent the minimum and maximum of the data. 
doi:1 0.1 371 /journal.pone.0080273.g003 
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Figure 4. Comparison of the distributions of task-induced modifications of global statistical network characteristics. Box-plots (see 
Figure 3 for details) of the relative deviation of the average shortest path length L* (left) and of the clustering coefficient C* (right) during the 
learning tasks (T[ and r') from the respective values during the baseline recording Tj, 
doi:10.1371/journal.pone?0080273.g004 



path length-of functional brain networks that we derived from 
ongoing multichannel EEG data recorded from 33 subjects during 
a neuropsychological test schedule. Despite the fact that learning 
and memory-related processes involve structures of the lateral and 
mesial temporal lobes, whose dynamics may not be directly 
accessible with non-invasive EEG recordings, we observed 
differential learning-related modifications of the networks' global 



statistical properties. While there were no detectable modifications 
of the average shortest path length, networks attained a 
significandy higher clustering coefficient during intentional learn- 
ing as compared to incidental learning and to the resting state. 
Interestingly, modifications during intentional learning even 
allowed us to predict the subsequent recall performance: the more 
pronounced the modifications of the clustering coefficient the 
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Figure 5. Relationships between task-induced modifications of global statistical network characteristics and recall performances. 

Scatterplots of relative deviations of the average shortest path length L* (top) and of the clustering coefficient C* (bottom) during incidental T[ (left) 
and intentional learning 7j (right) from the respective values during the baseline recording and subsequent recall performances N\ and N 2 . A 
significant correlation could only be observed between the relative clustering coefficient C* T , and the number N2 of recalled words during T' 2 (linear 

regression is represented with a solid black line). 
doi:1 0.1 371 /journal.pone.0080273.g005 
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higher the recall performance. An elevated clustering coefficient 
might be indicative of an increased occurrence of small but tighdy 
connected groups of network nodes. Identification of these groups 
(or communities) using appropriate analysis techniques [88,89] 
may aid in further exploration of the biological relevance of 
communities [29] and may provide novel insights into the 
relationship between the topology of functional brain networks 
and cognitive functions. 

Our findings are in line with some previous reports on task- 
related modifications of global properties of functional brain 
networks and their relationship to task performance (e.g. 
[19,24,29,35]) and support the view that higher cognitive 
functions-such as learning- and memory-related processes-result 
from dynamic interactions of distributed brain areas operating in 
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